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Recent  studies  of  two-cine- sc  ale  linear  systeas  have  used  a slaple  transformation  of  vari- 
ables to  block-diagonalize  the  systea  so  that  fast  and  slow  nodes  are  decoupled.  A matrix 
used  In  this  transformation  satisfies  an  algebraic  Rlccad  equation.  In  this  study  It  Is 
shown  that  this  transformation  is  unique  and  Its  elgenstructure  Is  described.  Also,  a new 
method  of  computing  the  decoupling  transformation  Is  presented  end  demouatreted  using  the 
16th-order  linear  model  of  a turbofan  engine. 

INTRODUCTION 

Consider  the  system  of  linear  first-order  differential  equations 

i - Ax  (1) 

where  x Is  an  n-dlmenslonal  vector  and  A Is  a real,  constant  matrix.  As  proposed  by 
Chow  and  Kokotovlc  (1976),  the  system  (1)  will  be  called  two-time- scale  If  the  eigenvalues 
of  A,  represented  as  A (A),  can  be  separated  by  absolute  value  Into  nonempty  sets  S and 
F so  that 

l»tl  « |fj!  for  all  st  6 s and  f^  « F.  (2) 

Let  n^  be  the  number  of  eigenvalues  In  S,  and  order  them  so  that 

|mtl  i l*t+1l . 1 • 1,  2,  ...  nl  - 1 

Ifjl  < l*j+1l.  i - 1.  2.  ...  n,  - 1, 


"2  * l’ 


WIIWIV  I»2  I*  U jS 

For  any  two- time-scale  system  the  ratio 


will  define  the  systea  small  parameter  needed  to  measure  the  system's  time- scale  separation. 
Such  syst«s  are  also  described  as  stiff,  ill-conditioned  or  singularly  perturbed. 

Given  the  system  with  n.  slow  modes  and  n,  fast  modes,  partition  the  state  vector  x 
into  subvectors  x.  and  x.  of  dimension  n.  and  n,,  and  partition  the  A matrix 
accordingly  so  (1)  can  be  wfltten  1 L 


ll 

AI2 

‘21 

A22 

The  following  two  step  transformation  has  been  used  by  Chang  (1972)  and  Kokotovlc  (1975)  to 
first  reduce  (4)  to  block- triangular  form  and  then  to  block-diagonal  form.  If  the  n,  * n 
matrix  L satisfies  the  nonsymmetrlc  algebraic  Rlccatl  equation  2 1 
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where  Q ■ K *.  The  H,  J end  Q matrices  are  perclclonad  compatibly  with  Che  parti- 
tioned system  (4). 

If  the  A aatrlx  la  nondefeccive,  l.e.(  haa  a full  aet  of  linearly  Independent  elgenvac- 
tora,  chan  J la  diagonal  with  elgenvaluea  appearing  along  Che  diagonal.  The  M aatrlx 
la  the  nodal  aatrlx  vhoaa  coluana  are  the  correapondlng  (right)  elganveccora  of  A.  If  A 
la  defective  then  J contalna  one  or  aore  Jordan  blocks  and  aoae  of  the  collates  of  M are 
generalised  eigenvectors.  Q Is  the  aatrlx  of  reciprocal  basis  vectors  whose  rows  are  the 
left  eigenvectors  and  left  generalised  eigenvectors  of  A.  The  existence  and  uniqueness  of 
M and  J are  discussed  In  many  linear  algebra  texts  such  as  Stewart  (1973) . 

For  the  case  where  A has  multiple  eigenvalues  It  can  be  difficult  to  get  a unique  repre- 
sentation of  M.  However,  If  1 la  an  eigenvalue  vlth  multiplicity  a,  there  la  a unique 
a-dlaenslonal  aubspace  spanned  by  the  eigenvectors  and  generalised  eigenvectors  correspond- 
ing to  X . 

Naraslahaaurthl  and  Wu  (1977)  and  Valjkn  (1977)  have  described  the  elgens  true  cure  of  the 
solution  to  the  algebraic  Rlccatl  equation  (3).  Those  results  are  given  below  as  Theory  1 
with  a simplified  proof.  Potter  (1966),  Martensson  (1971)  and  Rue  era  (1972)  presented 
slnllar  results  for  the  syametrlc  algebraic  Rlccatl  equation  which  arises  froa  the  linear 
regulator  problaa. 

In  Theorem  1 we  need  not  asaae  chat  A Is  two-tlme-ecale,  but  only  chat  partition  (16)  la 
possible.  That  Is,  that  Che  Jordan  aatrlx  J can  be  partitioned  Into  subaatrlces  J. , J, 
without  splitting  a Jordan  block.  Also,  the  L aatrlx  aay  be  complex.  1 2 

Theorem  1.  Matrix  L satisfies  the  algebraic  Xlccatl  equation  (3)  If  and  only  If  L 
satisfies  the  linear  equation 

“ll  * -"21  OM 

for  A - M J Q partitioned  as  (16). 


Proof:  Assume  first  that  (3)  la  satisfied.  Rewrite  (5)  as  L(A,,  - AjjL) 


and  let  (A^  - A^L)  have  Jordan  form  XGX*1.  Setting  T • -U,  It  follows  that 


11 


Poac-aulclply  by  X to  obtain 


"11 


*12L 


.-1 


A21  " *22L  ‘ "Ljrat 


*11*  + A12Y  " XG 
A2i*  + A22Y  - YG 


X X 

A - G. 

-YJ  LY. 

Thus  the  diagonal  elaments  of  the  Jordan  fora  C are  n^  of  Che  n eigenvalues  of  A,  and 

[:] 

are  a.  corresponding  eigenvector*  of  A.  Also  LX  - -Y,  completing  the  fir*t  half  of  the 
proof . 

Suppose  LMU  - -M21.  We  need  to  shov  that  the  lover  left  block  of  the  product 


"l  o' 

*11  *12 

' I o' 

*ir*i2L  *12 

L I 

. *21  *22  . 

-L  I 

U11+A21-LA12L-A22L  a22+LA12  ^ 

is  aero. 


Since  H is  full  rank,  the  matrices 
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Mu  and  Mu 


art  both  full  rank  and  L can  ba  rapraaancad  at  -M  M . Also,  since  Q M ■ I,  It  fol- 
lows that 


2111 

so  the  1 aatrlx  can  alao  ba  rapraaantad  as 


+ 922M21  * 0 


922L  ’ 921 


or  L-Q2‘Q21. 


(16) 


Writ*  A In  the  Jordan  form  and  complata  tha  product  so  that 


-M21Mli  1 
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M11J1MU  MUJl912+M12J2922  ■ 

"!1  *12  ' 

. 0 922J2922 
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complatlng  cha  proof. 


In  tha  final  expression  v*  uaa  the  idantltlaa 

w-1 


11 


911  " 912922921 


and 


922  ■ M22  - MaM'liM12- 


As  this  chao ram  illustrates,  if  L .satisfies  (15),  than  tha  eigenvalues  and  eigenvectors  of 
B.  and  B.  are  M..,  J.  and  Qli»  r«»P«ctiv*ly . Tha  casa  of  generalized  eigenvec- 

tors causes  no  difficulty1^  tha  proof  sine*  any  Jordan  block  appearing  In  G alao  appears 
In  J^,  and  vice  versa. 

Tha  conditions  which  Insure  that  L la  real  also  follow  quite  readily  from  this  theorem. 
That  Is,  If  A has  complex  eigenvalue  X appearing  In  J. , than  X and  X * must  only 
appear  In  J..  Under  this  restriction  there  will  always  exist  a nonsingular  complex  matrix 
C such  that1 


H 

X 
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Is  real.  If  L satisfies  LM. . - -M_.  Chen  L also  satisfies  LM'  • M’  so  that  L la 
real.  11  21  11  11 


The  general  form  of  the  block-diagonalization  transformation  (12)  can  be  determined  from 
Theorem  1.  Comparing  the  Jordan  forms  A • M J Q and 


UNIQUENESS  OF  THE  L AM)  K MATRICES 


I 


A*  shown  by  Naraslmfaamurthl  and  Wu  (1977),  only  undar  vary  raatrlctlva  conditions  will  tha 
algabralc  Rlccatl  aquation  hava  a unlqua  solution.  Namely,  only  whan  tha  quadratic  tan  la 
taro  and  tha  aquation  bacoaaa  llnaar.  Howavar,  thara  la  a unlqua  solution  which  decouples  a 
two- tins- ac ala  system  with  alow  and  fast  nodaa  aaparatad.  Thla  la  provan  In  tha  following 
theorem. 


Theory  2.  For  a glvan  two-tlaa-scala  ayatan  A thara  la  a unlqua  dscoupllng  matrix  L 
satisfying  both  (5)  and  (11). 

Proof:  As sums  that  L and  L'  both  satisfy  (5)  and  (11).  Applying  Thaoram  1 thara  axlat 
full  rank  matrlcaa 


' X ' 

* X' 

and 

. Y _ 

Li’J 

such  that  LX  • -Y  and  L'l'  * -Y ' . Moreover,  tha  coluams  ara  (ganaralizad)  alganvactora 
of  A corresponding  to  alganvaluas  S.  Sines  thara  la  a unlqua  n .-dimensional  subapaca 
spanned  by  thaaa  (generalized)  eigenvectors,  there  exists  matrix  such  that 


Post-multlply  L'X'  - -Y’  by  C Co  obtain  L'X  - -Y.  Aa  shown  In  tha  proof  of  Theorma  1, 
* exists.  Thus  L'  • L “ -YX  completing  tha  proof. 


Although  tha  decoupling  matrix  L la  unlqua  for  a glvan  two-tlae-scale  ayatan,  there  ara 
generally  many  real  solutions  to  tha  algabralc  Rlccatl  equation  (5).  Thus  any  effort  to 
obtain  L should  avoid  these  other  solutions.  In  a special  caaa  tha  nuaber  of  distinct 
solutions  can  ba  enumerated  by 


Corollary  1.  If  A has  distinct  real  alganvaluas,  than  thara  ara  ( ° I distinct  real 
solutions  of  (5)  corresponding  to  tha  number  of  distinct  ways  of  \ "l|  partitioning  tha 
n eigenvalues  of  A Into  sots  of  n^  and  n^  members. 


Proof:  Since  M la  real,  all  solutions  of 
l*  X 1 Prov«  chat  this  Implies 
Tha  linear  Independence  of  tha  eigenvectors 
valuaa  ‘(AjfA^)  u *(*11-*1!1''>  Implies 


(5)  ara  real.  We  will  assume  that  X(A,.-A,,L) 
1H'.  Let  LX  - -Y  and  L'X'  - -Y'11  12 

corresponding  to  tha  n,  + 1 or  more  elgen- 
that  tha  augmented  matrix 


X X’ 

. Y Y*  _ 


(18) 


has  rank  at  least  ^ + 1.  Than  L - -YX*1  and  L'  - -Y'  (X')*1.  If  we  assume  that  L - l' 
so  that  YX*1  - Y'(X')*1  and  Y - Y'U')*1!,  than 


’ X ' 

’ X 

■ 

’ X'  ‘ 

_ Y _ 

. Y’tt')*^. 

. Y'  - 

which  Implies  that  the  n,  columns  of  tyl  are  llnaar  combinations  of  tha  columns  of 
and  tha  rank  of  (18)  equals  n, . a contradiction.  Thus  L i L'  and  tha  proof  Is 
complete.  1 


X 

Y' 


Tha  uniqueness  of  tha  X matrix  la  Insured  by  tha  two-tlma-acala  property.  As  shown  by 
Cantmachar  (1959),  tha  llnaar  aquation  (8)  has  a unlqua  solution  provided  8,  and  B„  hava 
no  cnauu  alganvaluas.  1 ‘ 

THE  CHOICE  OF  SLOW  AM)  FAST  VARIABLES 

An  Important  question  which  arlaaa  In  tha  decoupling  and  reducar-ordar  modeling  of  large 
lin**r  systems  la  how  should  stats  variables  x ba  reordered  so  that  states  which  contain 
predominately  slow  soda  ara  placed  In  x . If  wa  raordar  tha  state  variables  than  tha  rows 
of  tha  modal  matrix  it  ara  reordered  inthe  same  manner.  From  our  previous  results  It  la 
apparent  that  *Ll  la  approximately  proportional  to  lM,.l  and  Inversely  proportional  to 
IK11  . So  If  wa  can  choose  an  ordering  which  minimizes  Im^I  » this  will  insure  that  tha 


i . 
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I 

[ < 


■low  aod*«  appsar  primarily  In  varlablaa  x. , and  ILl  will  ba  small.  In  a similar 
fashion,  Cha  X matrix  can  bs  written  as  1 


•*1*12 


so  chac  choosing  an  ordering  of  scats*  which  minimizes  IM^I  will  pise*  primarily  fasc 


varlablaa  in  x.  and  ylald  a small  IKI. 


This  procaaa  la  darlflad  by  introducing  alow  mode  and  faat  soda  coupling  ratios  p 

pr 


'M 


lH12l 


Hara 

antrlaa 


rafara  to  tha  Euclidean  norm  and,  whara  necessary,  taka  tha  modulus  of  complex 
Either  p or  p.  can  ba  minimized  by  searching  over  all  orderings  of  tha  row 
indices  {1}  of  the  modal  matrix  M ■ [m  ].  If  it  should  happen  that  all  rows  of 


11 


21 


have  equal  norm,  chan  cha  minimum  value  of 


will  ba  (n,/n. ) 


1/2 


On  cha  ocher  hand,  if 


x can  ba  ordered  so  thac  • 0,  chan  Pm*m  0 and  the  iloi  modes  are  naturally  de- 


coupled from  cha  x.  states;  For  this  ordering  of  x,  cha  A matrix  will  ba  block- 


' ' was- 

triangular.  Likewise,  if  M 


triangular. 


12 


Chan  pf  is  zero  and  cha  resulting  A will  ba  block- 


D spending  upon  cha  structure  of  M,  cha  ordering  which  minimizes  p may  or  may  not  also 
minimize  p..  Depending  upon  whether  we  are  primarily  Interested  ln*studylng  cha  slow  be- 
havior or  fisc  behavior  of  cha  system,  we  could  choose  cha  ordering  of  variables  which 
minimizes  alcher  p or  p^. 


If  both  o and  p-  are  small,  wa  could  describe  cha  system  as  being  "weakly  coupled." 
This  definition  of  weakly  coupled  systems  is  similar  Co  one  proposed  by  Milne  (1965).  Ho 
ever,  his  definition  also  required  thac  u be  small.  Even  though  a system  is  strongly 
cwo- time-scale,  it  need  noc  therefore  be  weakly  coupled. 


Provided  one  is  able  Co  compute  eigenvectors,  the  following  procedure  is  suggested  for  re- 
indexing x,  A and  M before  obtaining  cha  decoupling  crane forma cion. 


Algorithm  1. 

1.  Normalize  che  columns  of 
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21 


(19) 


so  thac 


£ 

i-1 


U1 


- 1. 
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2.  If  (19)  contains  one  or  more  complex  eigenvectors  v , 
replace  v^  with  Im(v^).  1 


replace  with  Ra(v^)  and 


Evaluate  Che  norm  of  each  of  che 


r.  - ( r1  .?/'* 

j-1 


4.  Order  che  set  (r.)  so  chat 
the  permutation  matrix 


> r 

*1“  °2 


rows  of  (19)  snd  1st 
1 - 1,  ...  n. 

».  r snd  prsaultlply  x,  A and  M by 


[S  •*.  - \] 


(20) 


wtnrs  a Is  ths  1-th  colissn  of  ths  n * n ldsntlty  nstrlx. 


9.  Poat-aultlply  A by  ths  Invars,  pamutatlon  nstrlx 


) 


where  [•]* 


Indicates  trsnsposs. 


This  gsnsrsl  method  for  systsmatically  picking  slow  variables  should  bo  vary  useful  in  appli- 
cations as  discussed  by  Tenekctsls  and  Sand ell  (1976).  In  general  not  all  eigenvectors  are 
required,  but  only  those  corresponding  to  the  slow  nodes.  However,  one  question  which  needs 
further  exploration  is  how  scaling  of  the  state  variables  influences  this  method. 


COMPUTING  THE  L AM)  K MATRICES 


Kokotovic  (1975)  proposed  a method  of  computing  the  L matrix  starting  with  initial 
approxlnatlon 

Lo  “ *22*21 

and  ualng  tha  iteration 

Ll+1  " *22<Ll(*U  ” *12L1)  + *21)- 


(21) 

(22) 


He  proved  that  If  a certain  criterion  Involving  matrix  norma  la  eatlafied  than  thla  ltera- 
tion  will  converge  to  L.  Aa  he  noted,  the  criterion  is  quite  restrictive  and  is  not  satis- 
fied by  many  typical  two- time-scale  systems. 


If  efficient  eigenanalysis  programs  are  available,  then  a computational  procedure  for  find- 
ing accurate  dinar ical  L and  iC  matrices  based  on  the  results  of  Section  2 is  quite 
straightforward.  Examples  of  such  programs  are  EISPACK  from  the  Argo one  National  Laboratory, 
described  by  Smith  at  al.  (1976),  or  EIGRP  from  the  International  Mathematics  and  Statistics 
Library  (IMSL),  Houston,  Texas.  These  eigenanalysis  programs  are  well  documented  and  easy 
to  implement  on  mainframe  computers.  Typical  execution  times  for  finding  all  eigenvalues 
and  eigenvectors  of  matrices  of  order  40  and  80  are  .66  and  4.6  seconds  using  an 
IK  370/195  computer,  cf.  Smith  et  al.  (1976).  Both  programs  are  based  upon  the  QR  algor- 
ithm, cf.  Wllkenaon  (1965). 


Sine,  Che  L matrix  can  b a up r.»«d  aa  .lth.r  -M-.M'J  or  , not  .11  eigenvector, 

of  A need  be  evaluated  to  solve  for  L.  Two  procedures  could  oe  used,  one  for  the  case 
°1  — • °2  and  *®°th*r  f°r  **1  > *»2*  *or  ni  — n2_  obtmin  ai  (tight)  eigenvectors  corre- 
sponding to  the  slow  eigenvalues  S and  solve  the  set  of  linear  equations  (15).  Such 
equations  can  be  efficiently  solved  by  some  Gaussian  elimination  scheme  such  aa  LU  decompo- 
sition, e.g.,  Stewart  (1973).  For  n.  > n_,  find  the  n?  left  eigenvectors  (Q?TQo7) 
corresponding  to  the  fast  eigenvalues  F Ind  solve  the  linear  system  (16).  In  the  event 
that  complex  eigenvectors  sra  encountered,  apply  step  2 of  Algorithm  1 to  eliminate  complex 
entries  in  (15)  or  (16). 


The  accuracy  of  the  matrix  L obtained  from  (15)  or  (16)  can  be  evaluated  by  substituting 
it  back  into  the  Riccatl  equation  (5)  and  evaluating  the  residual  error  matrix 


■ Lo*ll  + *21 


Lo*12Lo  ‘ *22Lo’ 


If  R la  Judged  co  b.  not  aufflclently  amall,  1 ....  IR  I > c for  some  null  number  c, 
than  the  accuracy  of  L can  ba  lnprovad  by  an  ltaratlva0 procedure.  If  L la  a reaaon- 
abla  approximation  to  L than  tha  matrix  (A„  + LA)  hae  large  elgenvSluea,  la  non- 
alngular,  and  augfeata  tha  Ilnur  ltaratlon  11  0 

Ll+1  * (*22  + Ll*12)  1(L1*11  + *21)- 


Thla  ltaratlon  can  ba  r .written  aa 
Algorithm  2. 

1.  Obtain  initial  approximation  Lo  from  wither  (IS)  or  (16)  and  aat  i • 0. 

2.  Evaluate  »t  • L^  + A^  - L1A12L1  - A^  and  atop  if  IR^  < t. 

3.  Solve  (A22  ♦ L1A12)D1  - Rt  for  Dj  and  let  Li+1  - Lt  + 

4.  Let  1 • i ♦ 1 end  go  to  2. 

It  esn  be  shown  thet  for  L close  to  L,  the  proposed  iteration  converges  with  epproxi- 
mete  rate  u.  That  is. 
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*m  • <L1  ♦ VAU  ♦ *21  - (ll  * Dl>*12ai  * V ' *22at  + Dl> 

Rl+1  " *1  + D1(A11  " A12L)  * (*22  * L1A12)D1  * °1A12D1 


and  caking  cha  corractlon  nacrlx  Dj  aa  In  acap  1 abova 

Ri+1  “ <A22  * L1A12)  Ri(All  ‘ A12Ll  * A12Di) 

• (A22  + “ *12^1+1^  * 

Tha  largaac  elganvaluea  of  (A22  ♦ LA.,)”1  and  (A^2  • AjjL)  ara  1/f  l 
if  L is  close  co  L(  Chon 


Rml  . N 

"PIT’TM 


u « l. 


a . ao  that 
"1 


The  convergence  of  che  iceracion  defined  by  Algorithm  2 could  be  Insured  by  Including  a line 
search  In  seep  3.  Thac  Is,  If  lR. .1  > IR.I,  Chen  lec  L.^.  - L ♦ hD  end  search  on  h 
undl  *R  I Is  minimized.  However,  c his 1 safeguard  procedure  was  noc  needed  In  Che 
example  which  follows. 

This  procedure  for  obcalnlng  numerical  solutions  Co  che  nonayametrlc  Rlccatl  equation  has 
some  similar  ides  co  one  proposed  by  Farrar  and  D1  Pietro  (1977)  for  solving  che  syMetrlc 
Rlccacl  equation  arising  In  che  linear  quadratic  regulator  problem  in  chat  both  procedures 
are  initialized  using  eigenvectors.  However,  cheir  Iterative  correction  lnvolvee  a bi- 
linear equation  whereas  Algorithm  2 Involves  a linear  Iterative  correction. 

Algorithm  3 listed  below  can  be  used  to  compute  the  R matrix.  Based  on  che  previous 
argument  that  the  eigenvalues  of  B*  are  all  larger  than  those  of  Bj,  It  can  be  shown 
that  the  rate  of  convergence  will  approximately  equal  che  system  small  parameter. 

Algorithm  3. 

1.  Sec  K • 0,  R ■ -A.,,  and  1*0. 

o o 12* 

2.  Solve  DjB^  ■ RA  for  and  let  ♦ D^. 

3.  Evaluate  R1+1  - - B1x1+1  ♦ mnd  *COp  ,R1+1*  — c# 

4.  Set  1 - 1 + 1 and  go  to  step  2. 

An  advantage  of  this  approach  is  that  B.  is  unchanged  in  the  Iteration  so  that  it  can  be 
decomposed  once  into  LU  form  and  then  used  repeatedly  In  step  2 to  find  D^. 

To  summarize,  the  proposed  procedure  for  computing  the  decoupling  transformation  for  a two- 
time-scale  system  Involves  che  following  steps. 

1.  Obtain  approximate  eigenvalues  X (A)  and  choose  the  number  of  slow  modes  n..  Note 
that  different  choices  of  n^  are  possible  to  suit  the  needs  of  the  study. 


2.  If  n-<  n.  obtain  approximate  eigenvectors 
vector.  [Q2i  Q22l - 

3.  (Optional)  Apply  Algorithm  1 to  placs  slow  variables  in  x^. 

4.  Eliminate  Imaginary  parts  as  needad,  and  solve  either  (IS)  or  (16)  for  L . 

o 

5.  (Optional)  Iaprovn  chn  accuracy  of  Lg  by  Algorithm  2. 

6.  Find  t by  Algorltla  3. 

Thin  procadur.  ha.  boon  auccaaaf ully  applied  to  flv,  nodala  of  phynlcal  .Titan,  with  ord.r. 
ranging  fro*  A to  32.  In  all  caaaa  AlgorlthBa  2 and  3 converged  with  approximate  rate 
u and  It  waa  verified  that  the  eigenvalue  decoupling  condltlone  (11)  ware  eatlaflad.  How- 
ever, for  tha  aake  of  brevity,  only  one  of  theae  axanplaa  will  be  preaented  here. 


11 

"21 


Otherwise,  obtain  left  eigen- 


TURBOFAN  EXAMPLE 


This  example  Is  s 16th  order  mods!  of  s curbofsn  engine  which  was  the  theme  problem  for  s 
recent  conference  on  control  of  linear  multivariable  systems,  i.e.,  Sain  (1977).  This  model 
is  the  linearization  of  a detailed  nonlinear  simulation  of  the  turbofan  at  the  sea  level 
static  Intermediate  thrust  operating  point.  Two  of  the  state  variables  are  shaft  speeds, 
three  are  pressures  and  eleven  are  temperatures.  The  eigenvalues  for  the  linear  system  are 
-0.65,  -1.90,  -2.62,  -6.72±jl.31,  -17.8±j4.8,  -18.6,  -2l.3±j0.8,  -38.7,  -47.1, 

-50.7,  -59.2,  -175.7,  and  -577.0.  For  n.  ■ 15,  5 and  3,  the  resulting  small 
parameters  are  0.304,  0.371,  and  0.383.  1Since  these  values  are  not  particularly  small 
relative  to  one,  we  might  call  this  system  weakly  two-tlme-scale. 

The  computations  were  performed  at  the  University  of  Arizona  Computing  Center  on  both  the 
DEC  10  and  CYBER  175  computers  to  compare  the  effects  of  computer  word  length  and  compare 
execution  times  for  time  sharing  against  batch  processing. 

The  performance  of  Algorithms  2 and  3 for  the  two  cases  n*  • 3 and  n.  • 5 is  illustrated 
in  Figures  1 through  4.  The  state  variables  were  not  reordered  for  this  computation,  i.e.. 
Algorithm  1 was  not  employed.  As  shown,  in  all  cases  the  convergence  rate  approximately 
equals  the  small  parameter.  The  corresponding  computation  times  are  listed  in  Table  1.  The 
six  orders-of-magnltude  difference  between  the  minimum  residual  errors  obtslned  with  the 
DEC  10  and  CYBER  computers  can  be  directly  attributed  to  their  36  and  64  bit  word  lengths. 
Both  EISPACK  and  IMSL  were  tested  in  the  computation  of  LQ.  Their  performance  in  terms  of 
eigenvalue  accuracy  was  found  to  be  very  similar,  with  1SML  yielding  somewhat  shorter  execu- 
tion times.  The  execution  times  In  Table  1 were  obtained  using  EISPACK. 

T 

The  effects  of  reordering  the  state  variables  before  partitioning  x into  (x. ,x.)  are 
illustrated  in  Figure  5 for  the  case  n.  • 5.  As  shown,  the  choice  of  ordering  does  not  in 
this  case  affect  the  rate  of  convergence  of  Algorithm  2,  but  strongly  influences  p and 

ILI.  * 

The  computation  of  the  L matrix  using  alternate  initial  values  L is  described  in  Figure 
6.  It's  interesting  to  note  that  Algorithm  2 succeeds  in  converging  to  the  correct  solution 
when  initialized  with  L • A , i.e.,  equation  (21),  or  with  L ■ 0.  This  demon- 
strates that  in  some  casSs  tlM  decoupling  transformation  can  be  obtained  without  the  use  of 
eigenvectors.  However,  the  L matrix  obtained  using  eigenvectors  provides  a much  more 
reliable  initial  value  for  Alg8rlthm  2. 

CONCLUSIONS 

This  work  considers  the  properties  of  a transformation  used  by  Chang  (1972)  and  Kokotovlc 
(1975)  to  transform  linear  systems  into  block-diagonal  form  with  slow  and  fast  modes  de- 
coupled. The  elgenstructure  of  this  transformation  is  described,  and  it  is  shown  that  for 
any  given  two-tlme-scale  linear  system  the  decoupling  transformation  is  unique. 

Ceneral  conditions  for  the  existence  of  solutions  of  the  nonsymmetrlc  algebraic  Rlccatl 
equation  are  established  which  are  more  general  than  existence  conditions  previously  stated 
by  Naras lmhamur thl  and  Wu  (1977).  It  Is  also  shown  that  when  the  corresponding  A matrix 
has  distinct  real  eigenvalues,  the  distinct  real  solutions  to  the  Rlccatl  equation  can  be 
specifically  enumerated. 

A new  method  of  computing  the  decoupling  transformation  is  proposed.  With  this  method  the 
number  of  slow  modes  can  be  chosen  to  meet  the  needs  of  the  study  provided  multiple  or  com- 
plex eigenvalues  are  not  separated.  The  method  includes  algorithms  for  selecting  slow  sys- 
tem variables,  and  computing  accurate  L and  K matrices.  Although  system  eigenvectors 
are  used  in  this  method,  not  all  eigenvectors  need  be  found,  and  highly  accurate  eigenvec- 
tors are  not  required.  A 16th  order  linear  model  of  a turbofan  is  decoupled  by  this  method 
to  demonstrate  computation  times,  convergence  rates,  and  algorithm  stability. 
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Figure  1.  Computation  of  L for  the 

Figure  2.  Computation  of  K for  the 

case  - 3,  u - 0.383. 
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Figure  3.  Computation  of  L for  the 
case  n^  ■ 5,  u • 0.371. 
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Figure  4.  Computation  of  K for  the 
case  n1  - 5,  u - 0.371. 
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Figure  3.  Computation  of  L with  reordared  Figure  6.  Computation  of  L with 
state  variables  (nx  - 5).  | alternate  forma  of  L_ 


Xj  Indices  I Li  o* 

-(1,2, 3, A, 5)  3*2.  .260 

-(1,2,5,6,10)  213.  . 2U8 

-(1,2,9,11,14)  39.  .187 

^chosen  by  Algorithm  1 


(nr  - 5). 
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Table  1.  Typical  Execution  Times 
(CPU  seconds) 


Time  to  compute  L 


Case  n^  • 5 

CYBER  ] DEC  10 


.109 


Case  n^  • 3 


CYBER  DEC  10 


.599 


Time  per  iteration 
(Algorithm  2) 


Time  per  iteration 

(Algorithm  3) 


.015  .126 


.021  .171 


.012  .091 
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